Experimental study of isovector spin sum rules 
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We present the Bjorken integral extracted from Jefferson Lab experiment EGlb for 0.05 < < 
2.92 GeV^. The integral is fit to extract the twist-4 element /f " which appears to be relatively 
large and negative. Systematic studies of this higher twist analysis establish its legitimacy at 
around 1 GeV^. We also performed an isospin decomposition of the generalized forward spin 
polarizability 70. Although its isovector part provides a reliable test of the calculation techniques 
of Chiral Perturbation Theory, our data disagree with the calculations. 

PACS numbers: 13.60.Hb, 11.55.Hx,25.30.Rw, 12.38.Qk, 24.70.+S 



INTRODUCTION 

The Bjorken sum rule fi\ relates an integral over the 
spin distributions of quarks inside the nucleon to its axial 
charge. This relation has been essential for understand- 
ing the nucleon spin structure and establishing, via its 
Q^-dependence, that Quantum Chromodynamics (QCD) 
describes the strong force when spin is included. The 
Bjorken integral has been measured in polarized deep 
inelastic lepton scattering (DIS) at SLAG, GERN and 
DESY [3|-[3 ^'id at moderate four- momentum transfer 
squared at Jefferson Lab (JLab) Q, see e.g. Ref. [§| 
for a review. The variable is inversely related to the 
space-time scale at which the nucleon is probed. In the 
perturbative QGD (pQGD) domain (high Q^) the sum 
rule reads [iQl] : 



rr"(Q2) ^ r dx {g{{x, Q') - g^x, Q^)) = 
Jo 
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where g\ and g" are the spin-dependent proton and neu- 
tron structure functions, gA is the nucleon axial charge 
that controls the strength of neutron /3-decay, as{Q'^) 
is the strong couphng strength and x = Q'^ /2Mv, with 
V the energy transfer and M the nucleon mass. The 
bracket term (yLt2, known as the leading twist term) is 
mildly dependent on due to pQGD soft gluon radi- 
ation. The other term contains non-perturbative power 
corrections (higher twists). These are quark and gluon 
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correlations that need to be understood to describe the 
nucleon structure away from the large limit. The 
Q^-dependence of //2i(Q^) is calculable in principle from 
pQGD. In practice, this has been done for ^2 and ^4 
only [llj. We stress that, as is almost always the case 
with pQGD, although the Q^-dependences are known, 
the absolute values of ^2 and ^4 are unknown and need 
to be measured or computed by non-perturbative means. 
Besides its contribution to establishing pQGD (at high 
Q^), the Bjorken sum rule can be used to extract higher 
twists, to check lattice QGD calculations (at moderate 
Q^), and to test effective theories of the strong force (at 
low Q^). In addition, Bjorken sum data and phenomeno- 
logical models at lower can be described with a nearly 
constant "effective strong coupling" as,gi [13, [13. The 
lack of Q^-dependence of Us^gi opens new avenues for 
non-perturbative QGD calculations using the AdS/GFT 
correspondence jl4| . 

The elastic contribution to the Bjorken sum is usually 
not included because the generalized Bjorken sum rule is 
derived at large where such contribution is negligible. 
Furthermore, the Bjorken sum rule naturally connects 
to the Gerasimov-Drell-Hearn sum rule [isj in which the 
elastic is inexistent. Gonsequently, when presenting the 
experimental measurement of the Bjorken sum, the elas- 
tic contribution will not be included. We refer to ref. 
for a discussion on whether to include or not the elas- 
tic contribution to the GDH sum rule. However, for 
higher twist analysis, all reactions should be included for 
a meaningful higher twist extraction [3]-[l3. Therefore, 
in the part of the paper discussing higher twist extrac- 
tion, the elastic contribution to the Bjorken sum will be 
added. 

In this paper, new data from the JLab GLAS EGlb 
experiment [20|- 22| taken on polarized proton and 



deuteron targets are used to extract the Bjorken inte- 
gral over an extended range: 0.05 < < 2.92 GeV^ 
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compared to the previous JLab range 0.15 < < 1.5 
GeV2 [a]. 



The extension down to = 0.05 GeV^ allows us to 
compare to Chiral Perturbation Theory (xPT) calcula- 
tions in a domain where the chiral approximation should 
be valid. The moderate range had been precisely 
measured The new data set, of equivalent precision, 
provides a useful check. In particular, it verifies the neu- 
tron results, which come mostly from ^He in Ref. Q and 
from the deuteron in this paper. At larger {> 1 GeV^, 
where Eq.[l]holds), higher twists can now be studied with 
a statistical precision typically improved by a factor of 2. 
Previous work 0] has shown the necessity of precise 
mapping at moderate (> 1 GeV^) because of the sur- 
prisingly small size of the overall higher twist effect. One 
might be tempted to lower the values at which the 
analysis is done (see Eq. [1]) but this is not reliable due to 
the fast 1/Q^'~^ rise of twist i contributions at low 
and to the increasing uncertainty of the evolution of the 
twist-2 parts. The main contributor at low to this 
uncertainty is the strong coupling constant as{Q'^)- 
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TABLE I: The measured (r^';;;;;^^ and total {V\~"t) Bjorken 
integrals for different points (in GeV'^). The experimen- 
tal systematic uncertainty cr^eas is given in the 3'^'* column. 
Total systematics uncertainty, including the low and large- a; 
extrapolations, (cr"^"*) and statistical uncertainty (ct**"*) on 
^I'tot given in the 5"* and 6**" columns. 



The Bjorken integral is advantageous compared to the 
individual moments V\ and F" because of simplifications 
arising from its non-singlct (p-n) nature: at moderate 
lattice QCD calculations are easier and more reliable 
because disconnected diagrams, which cannot be easily 
computed on the lattice, cancel out. At higher Q^, the 
(p-n) simplification provides a sum rule (the Bjorken sum 
rule) based on more solid grounds than the sum rules 
for individual nucleons (the Ellis Jaffe sum rules [l^ 
that necessitate additional assumptions). At low , 
the (p-n) subtraction cancels the A1232 resonance con- 
tribution which makes the xPT calculations significantly 
more reliable [35^. By a similar argument, the transverse- 
longitudinal polarizability 5lt [SI, a higher moment of 
spin structure functions, also provides a reliable test of 
xPT computations. (In that case, the A1232 contribu- 
tion is suppressed at low because the N-A transition 
is mostly transverse, making the longitudinal-transverse 
(LT) interference term very small.) Nevertheless, calcula- 
tions based on xP'P and data for 5lt on the neutron [i^l 
strongly disagree. This calls for more low studies, 
especially the yet unmeasured [ll]. The data dis- 
cussed in this paper were taken with a longitudinally po- 
larized target and hence cannot be used to extract S^j^j.. 
However, the generalized forward spin polarizability 70 
can be obtained and, just like the Bjorken integral, its 
isovector part 7g — 7g offers the same advantages as 5lt 
for checking the calculation techniques of xPT- We will 
also report on these results. 



BJORKEN SUM EXTRACTION AND 
COMPARISON WITH CHIRAL PERTURBATION 
CALCULATIONS 

The measurements of structure functions g\ and gf 



22l |. The data cover an in- 
3 GeV for 0.054 < < 



are described in Refs. [2C 
variant mass range up to W 
2.92 GeV^. Since experimental moments are integrated 
over a finite W range, the data have to be supplemented 
by models for large W . We used the model described 
in Ref. down to x = 0.001. This part is known 
from DIS experiments. The rest is determined using a 
Regge parametrization l8| which was compared to that 
of Bass and Brisudova [28[ and found consistent with it. 
A parameterization was also used to estimate the con- 
tributions between pion threshold (1.08 GeV) and 1.15 
GcV 

The Bjorken integral is obtained from V\ and as- 
suming: 

Fr" = 2F?-Ff/(l-1.5c.,), 

with the deuteron D-state probability ojd — 0.05 ± 
0.01 '2§|. The data are given in Table I (a more detailed 
table is given in [sO]) and shown in Fig. [TJ The elas- 
tic contribution (x = 1) is excluded. Data from SLAG 
E143 i\, HERMES f?], JLab GLAS EGla (proton and 
deuteron), and JLab Hall A E94010 (neutron from ^He) 
combined with GLAS EGla (proton) [§| are also shown 
for comparison. 

There is excellent agreement between the Bjorken inte- 
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FIG. 1: (Color online) The Bjorken integral r^'"((32). The 
solid blue circles give the results from this work with the 
horizontal band giving the systematic uncertainties. Other 
symbols show the data from experiments E143 0| (open dia- 
monds), E155 H (open star), HERMES [3 (open triangles) 
and JLab 0] (open circles and open squares). For those, the 
error bars represent the quadratic sum of the statistic and sys- 
tematic uncertainties. The gray band represents the leading- 
twist NNLO pQCD calculation. The curves correspond to 
XPT calculations ( [3^ . [sj) and phenomenological models 

( MM)- 



gral with the neutron extracted from the deuteron (filled 
circles and open squares) and from ■^He (open circles). 
The neutron spin structure functions extracted from the 
deuteron and from ^He agree at moderate and large Q"^ . 
However, for below a few tenths of a GeV^, nuclear 
effects beyond those accounted for in the convolution 
method employed to extract the neutron |23| may be- 
come large 2J] . Therefore, at low one needs both the 
deuteron and ^He data to ensure a reliable neutron ex- 
traction. Nuclear effects in the deuteron are weaker, but 
there is an unsuppressed contribution from the proton. 
On the other hand, ^He is more tightly bound, but the 
polarized proton contribution is largely suppressed. Con- 
sequently, the uncertainty due to nuclear effects is mostly 
of different origin in the deuteron and '^He, which makes 
the two nuclei complementary. The agreement between 
the deuteron and '^He results is also encouraging for the 
interpretation of the low '^He and the deuteron data 
(Q^ > 0.015 GeV^) that wiU be available shortly respec- 
tively from Jefferson Lab's Hall A [2^ and B j2i]. The 
data also agree well with the SLAG and HERMES exper- 
iments and with the two phenomenological models shown 
in Fig. [T] The model of Burkert and loffe [3l[ (contin- 



uous black curve) is a meson-dominance-based extrapo- 
lation of DIS data supplemented by a parametrization 
of the resonance contribution. The other model (Soffcr- 
Teryaev |32'|, dashed red curve) uses the smoothness of 
gi + 52 with to extrapolate DIS data at lower . 

At moderate , we observe a strong variation of the 
Bjorken integral, in contrast to the high region. To- 
gether with our data at the lowest points, the kine- 
matic constraint Fi ^ when ^ suggests a small 
Q^-dependence of F^~" at low as well. This would 
agree with the fact that the Fi slope at ~ is given 
by the generalized GDH sum rule which predict a small 
Q^-dependence. 

At low the data are consistent up to 0.2 GeV^ 
with the calculations of Bernard et al. [s^ and up 

to Q^~0.35 GeV^ for those of Ji et al. done in the heavy 



baryon approximation 3J| . The range of validity of the 
xPT calculations seems larger than of individual nucle- 
ons 0, [i^l possibly because the A1232 resonance is sup- 
pressed in the Bjorken integral [3^. This result, however, 
is not trivial: Good agreement was expected between 6lt 
and xPT results since the A 1232 is strongly suppressed 
at low Of for Slt- However, its measurement for the 
neutron 



27[ disagrees strongly with xPT calculations. 



To quantitatively compare with xPT calculations, we 
fit our results up to a maximum ranging from 0.30 
to 0.50 GeV^ (fits on lower ranges are imprecise, and 
higher data may lie out of the region of validity for 
xPT). We included the data from Ref. [sl in the fit. Our 
fit form is: 



8M 



(2) 



in which k is the anomalous moment of the nucleon and 
a and b are fit parameters. The first term in Eq. [2] stems 
from the Gerasimov-Drell-Hearn sum rule 0. We find 
a = 0.80 ± 0.07{stat) ± 0.23(s?;si) and b = -1.13 ± 
0.16(stoi)±0.39(syst) with x^/dof = 1.50. The term 
agrees well with the results from Ji et al. {a ~ 0.74) but 
not with those of Bernard et al. (a = 2.4). The fit un- 
derscores the importance of the term (not calculated 
yet in xPT). This was also noticed for Ff and Ff [22]. 

At high Q^, the leading twist pQGD calculation is 
given by the bracket term of Eq. [T] and is represented by 
the gray band in Fig.[T] It agrees reasonably well with the 
data. This implies that the total higher twist contribu- 
tion is relatively small even down to ~ 1 GeV^ where 
one would expect higher twist contributions to be signif- 
icant. Higher twists, which measure parton correlations, 
are weighted by l/Q^'"^^ (with t being the twist number) 
and are related to the confinement mechanisms and to 
scattering off coherent quarks. Because of these reasons, 
it was initially expected that higher twists would play an 
important role at < 1 GeV^. Higher twists can be 
positive or negative but there is no fundamental reason 
to expect a well-tuned cancellation of different terms in 
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the higher twist series that would make the overall higher 
twist contribution small. However, this seems to be the 
case experimentally, at least around 1 GeV^. One 

of the aims of the higher twist analysis reported here is 
to establish whether higher twists are intrinsically small, 
or whether the terms in the higher twist series conspire 
to cancel. 



HIGHER TWIST ANALYSIS 



The first higher twist correction term in Eq. [T] is 11| : 
~9~ 



p—n 



4dr" 



4/r") 



(3) 



where a2 and ^2 are known. They are given by moments 
of the leading twist part of gi and the twists 2 and 3 parts 
of g2- a2 = /q dx {x gi) and d2 = /g dx x^ {2gi + 3.92). 
The twist-4 term that we wish to extract is /!"""• 

To perform a higher twist analysis, the elastic contribu- 
tion (x — 1) to r^~" is added. The moment F^"" which 
includes the elastic contribution estimated from form fac- 
tor parameterizations [s^ is shown in Fig. [2l In Eq. [U 
as is computed up to next to leading order. A fit of polar- 
0.3 f 
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FIG. 2; (Color online) World data on the Bjorken integral, in- 
cluding the elastic contribution. The error bars represent the 
quadratic sum of the statistic and point-to-point-uncorrelated 
systematic uncertainties for the JLab data, and the quadratic 
sum of the statistic and full systematic uncertainties for the 
rest of the data. The continuous line is our three parameter 
fit in the g^-range from 0.66 to 10 GeV^. 



ized quark distributions [ll] yields " — 0.031 ± 0.010 
at g2= 1 GcV^, whereas d^'"" = -0.007 ± 0.010 is ob- 
tained from Ref. [3 and Ref. ^ evolved to 1 GeV^. The 
EGlb data on F^~", together with the world's data, can 
then be fit to extract /f " using Eqs.[T]and[31 To account 
for twists greater than rank 4, we include a coefficient 



/j.g~"'/(5^. For consistency, former data on F^~" were re- 
analyzed using the same model as used in this paper to 
extrapolate to low x. For both JLab data sets (Ref. Q 
and the present data), the point-to-point correlated un- 
certainties have been separated from the uncorrelated 
ones. The latter are added in quadrature to the statistical 
uncertainties. The correlated systematics are propagated 
independently, as is the uncertainty arising from as ■ The 
result of the fit done in the Q^-range from 0.66 to 10.0 
GeV2 is /f""(g2 = 1 GeV2) = -0.101 ± 0.027±0:!]f3 
with ne/M^ = 0.084 ± O.Ollibg:^;^. The first uncertainty 
is the quadratic sum of the statistical and the point-to- 
point uncorrelated uncertainties. The second one is the 
point to point correlated uncertainty. Comparing the val- 



ues of ",02 " and d^ " at 1 GeV^, we see that 
A^r" ~ 0-4/2^"" GeV2. The resuh for /f"" is plotted in 
Fig. [3] (square) along with the result from Ref. [H (tri- 
angle) and theoretical predictions (In addition to /2 and 
fie, the third fit parameter mentioned in Figs. [2] and [3] 
is ga, which was free to vary within its experimental un- 
certainty). As discussed in the introduction, only the 
Q^-dependence of /2 is known from pQCD. The abso- 
lute value can be computed solely from non-perturbative 
means and is difficult to obtain with Lattice QCD. For 
these reasons, only phenomenological models are avail- 
able for comparison with our results. 
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FIG. 3: fr"iQ = 1 GeV^) for the fits performed over the 
0.66 < < 10 GeV^ range for this study and Ref. 8] (JLab 
A/B). Calculations 40]- 44] are shown by the bands. Sum 
rule (1) refers to Ref 41]" and (2) to Ref g^. 

At = 1 GeV^, the leading twist term /i2~" and 
higher twist terms fJ-^~"' and /if ~" are of similar sizes but 
with alternating signs and with fJ^ " and /ig~" mostly 
canceling each other. 

To study the systematics associated with this higher 
twist analysis and to check the legitimacy of our proce- 
dure at low Q^, we conducted several tests: 

1. We repeated the fit for several ranges; 

2. We reiterated this work adding a lig~" /Q^ term 
to study the convergence of the twist series (the 
resulting /f " is shown in Fig.[3]by the solid circle); 

3. We investigated the dependence on the low x ex- 
trapolation using different Regge-based parameter- 
izations; 
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4. We extensively studied the stability of the fit for 
different choices of number of parameters and of 
ranges by using different models that reproduce the 
data reasonably well. We used ranges from 0.47 < 
0^ < 10 to 3 < Q2 < 10 GeV^ and we fit with 
functional forms with highest term from /^e/Q^ to 

All observations supports the validity of our extractions. 
See Ref. [s^] for details. 



derivation akin to that of the GDH sum rule: 

^"^^^io " (,-9i-^-52jdx, (5) 

where is the second spin structure function and a is 
the fine structure constant. Similar relations define the 
generalized longitudinal-transverse polarizability '■ 

ne\9LT(y.Q')~9tT(yM^)\^ (6) 



COLOR POLARIZABILITIES 

Combination of higher twist coefficients can be inter- 
preted in terms of color polarizabilities, which describe 
the response of the color magnetic and electric fields to 
the spin of the nucleon. The color electric and mag- 
netic polarizabilities (i^l, [3] are xe = § (2^2 + /2) 
and Xb = \ — j-i)- Using the value of /f"" ex- 
tracted from the fit with (5^jj„ — 0.66, we obtain Xe"" — 
-0.077 ± 0.050 and Xs"" = 0.024 ± 0.028. The point- 
to-point correlated and uncorrelated uncertainties on /2 
were added in quadrature. Our higher twist analysis 
yields I/2 "I ^ Mf""! feature predicted by mod- 
els nil and [H). Consequently x^" - f/f'" and 



3LT 



1 j-p-n 



ELECTROMAGNETIC POLARIZABILITY 

We now turn to the generalized forward spin polariz- 
ability 70. Spin polarizabilities characterize the coherent 
response of the nucleon to photons. They are defined 
using low-energy theorems in the form of a series expan- 
sion in the photon energy. The first term of the series 
comes from the spatial distribution of charge and current 
(form factors) while the second term results from the de- 
formation of these distributions induced by the photon 
(polarizabilities). Hence, polarizabilities are as impor- 
tant as form factors in understanding coherent nucleon 
structure. Generalized spin polarizabilities describe the 
response to virtual photons. The low energy theorem 
defining the generalized forward spin polarizability is: 



In 

{—)IrT{Q')>^ + lo{Q'y + Oii.''), 



(4) 



where qtt is the spin-flip doubly-virtual Compton scat- 
tering amplitude, and Itt is the coefficient of the 0{iy) 
term of the Compton amplitude which can be used to 
generalize the Gerasimov-Drell-Hearn (GDH) sum rule 
to non-zero We have Itt{Q^ = 0) = k/4. In 

practice 70 can be obtained from a sum rule which has a 



(51 + 92) dx. 



(7) 



where (?lt is the longitudinal-transverse interference am- 
plitude, and Ilt is the coefficient of the 0{v) term of 
the Compton amplitude. Details on the derivation of 
Eqs. lillTl can be found in [§] and The isovector 

quantity 7^— 7o eliminates the A 1232 resonance con- 
tribution [351, and therefore offers the same advantage 
as when comparing to calculations based on xPT- 
Higher moments are advantageous because they are es- 
sentially free of the uncertainty associated with the low 
X extrapolation. An isospin separation of 5lt or 70 may 
help us to understand why the xPT calculations fail to 
describe them. For example, the ^-channel exchange of 
axial- vector mesons (short range interactions), which are 
not included in the calculations, could be identified if one 
of the isospin components agrees with the xPP calcula- 
tions while the other disagrees. 

We formed 7q — 7q using the proton data from 
EGlb |22] and the neutron data from JLab experiment 
E94010 [27|. The ^He data [271 are more precise than 
the deuteron data [2^ that contain contributions from 
quasi-elastic and two-body break-up, which are not re- 
solved by the CLAS spectrometer but are large at low 
. (This difficulty prevented 7g from being obtained 
from the EGlb data l^j). EGlb goes to lower than 
E94010, but the coverage of E94010 is sufficient for our 
investigation. The resulting 7q — 7g is shown in Fig. [4] 
(top plot) together with the predictions from Bernard et 
al. at 0(P4) dl and Kao et al. &t 0{P^) Ex- 
perimental values are given in Table II. We also plot the 
result from the 2003 MAID model Ei. As is true for 
2^ 



and 7o [27[, xPT calculations disagree with 79"" 
as well. Clearly, the discrepancy seen for 7g and 79 can- 
not solely be due to the A1232 resonance. The MAID 
model, which provides a relatively good description of 79 
and 7q , disagrees mildly for their difference at the lowest 
point. Complementary to this study, we formed the 
isoscalar part 7o -l-7o and compared it to the data (Fig. [4] 
bottom plot). The gray band on the Bernard et al. result 
is due to the uncertainty from the A1232 resonance. The 
MAID model provides a good description, whereas the 
xPT-based calculations still disagree. A disagreement in 
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0.0040 
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0.0037 


0.0090 


0.9 


0.0217 


-0.0103 


0.0016 


0.0040 



TABLE 11: Isovector and isoscalar parts of the generalized 
forward spin polarizability 70. 



the xPT calculation of one of the isospin components of 
7o along with agreement for the other component might 
have allowed us to identify a missing piece, such as for ex- 
ample a short range interaction due to heavy mesons, in 
the xP'P calculations. However, the discrepancy between 
data and xPT calculations for both isospin components 
does not allow us to draw such conclusion. This suggests 
that the non-resonant background is responsible. 
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FIG. 4: The isovector 7q — 70 (top) and isoscalar 7q -|- 70 
(bottom) generahzed forward spin polarizabilities together 
with xPT-hased calculations and the results from the MAID 
model. The proton and neutron data are respectively from 
CLAS [2i and Hall A El 



SUMMARY AND CONCLUSION 

The Bjorken integral was extracted from polarized pro- 
ton and deuteron data for 0.054 < < 2.92 GeV^. The 
results for intermediate (the parton to hadron transi- 
tion domain) are consistent with previous JLab data in 
which the neutron information was extracted from po- 
larized '^He. This region exhibits a strong Q^-behavior, 
both from pQCD evolution and from some higher-twist 
effects. On the other hand, in the high-Q^ domain the 



Bjorken integral is rather flat. The data together with 
kinematic constraints at ^ also suggest a small 
Q^-dependence, in qualitative agreement with the gener- 
alized GDH sum predictions. 

At the lowest accessed by our data, xPP calcula- 
tions agree better with the Bjorken integral (an isovec- 
tor quantity in which the A1232 resonance does not con- 
tribute) than with moments on individual nucleons. This 
is not trivial since the xPP calculations fail to describe 
the generalized spin polarizability S]^x in which the A1232 
is also suppressed. 

Data on the generalized forward spin polarizability 
are not reproduced by the x^T^-based calculations 
even though the A1232 does not contribute. 

It is clear from previously published data on Slt and 
our analysis of 70 that the A1232 resonance contribution 
is not responsible for the discrepancy between data and 
calculations. The discrepancy between the xPT calcula- 
tions and the data occurs in all isospin channels, which 
makes it less likely that it is due to the contribution from 
heavier mesons in the chiral expansion. 

The low xPP regime has been recently mapped 
by two additional dedicated experiments in CLAS using 
polarized proton j4§| and deuteron targets [2^ and one 
in Hall A using polarized '^He 25|. These experiments 



p 
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will provide further precision tests of xPP calculation 
techniques. 

The moderate data (1 to 3 GeV^) allow us to ex- 
tract higher twist contributions and color polarizabil- 
ities. The twist-4 coefficient was found to be large: 
/f ~ -0.1 at Q2 = 1 GeV2 (compare to rf"" = 0.125, 
flP-" ^ 0.031 and d^~" = -0.007). The uncertainty 
on /f"" remains relatively large (w 70%); however, we 
have completed several systematic studies both with the 
existing data as well as simulated data (with no statis- 
tic fluctuations) that indicate our result is stable. The 
sign and magnitude of /I'"" agree with a recent analy- 
sis performed on gi directly [50|. The observation that 
higher twist effects on r^~" are small overall does not 
imply that the net higher twist effect on the structure 
function 5^"" is small at any x. It is important to 
study the x-dependence of the higher twists, as is done 
in Ref. (H^]. That I/2I is significantly larger than d2, and 
that /2 < 0, agrees well with the prediction of the two- 
scale model 4J| . Overall the net effect of higher twists is 
small, because of a cancellation between the twist 4 and 
twist 6 terms that are of similar sizes but opposite signs. 
This trend has also been seen for higher twist analyses 
done on the unpolarized structure function F2 [sH . This 
can be interpreted within a vector dominance framework: 
the oscillating signs arise from the development in series 
of the vector meson propagator oc 1/{Q^ — M"^) where 
Mm is the meson mass. 
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